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Loolution  of  the  Soimnerf eld-Gexl  Characteristic 
Equation  for  the  Stability  of  Couette  Plov/. 


1 t )  Introduction 


One  of  the  basic  problems  of  hydrodynamics  is  the  stabl 
llty  of  Couette  flow.  This  Is  shown  by  the  Importance  which 
since  ever  was  attributed  to  Couettes  observation  of  transi¬ 
tion  to  turbulence  [ 1  ]  as  well  as  to  a. I.  Taylors  investl- 
tiations  which  discovered  'the  occur^jncc  of  cellular  laminar 
vortices  at  certain  Ueynolds  numbers  |^2  j  .  The  intrinsic 
difference  of  these  two  observations,  once  transition  to  a 
random  eddy  motion  and  on  the  other  hand  tx'ansitlon  to  a 
laminar  vortex  pattern  is  illustrated  by  tho  two  limitiuj:, 
cases  in  which  they  occur,  the  transition  to  turbulence 
if  the  inner  cylinder  is  at  rest  and  the  outer  rotating: 
and  the  ijeneration  of  laminar  celluj.ar  vortices  in  the  re¬ 
verse  case  when  the  inner  cylinder  is  rotucii.j  and  the  outer 
at  rest.  Correspondiniily  the  mathejatical  treatment  is 
different.  V/hilst  for  the  rotatinis  outer  cylinder  Kayleitdis 
perturbations  [  3]  were  introduced  vxhich  still  are  the  oasis 
for  stability  investitiatlons  in  relation  to  the  generation 
of  turbulence,  for  the  rotatlnis  inner  cylinder  a. I.  Taylor 
introduced  perturbations  foxnning  rotational  symmetric 
cellular  vortices. 

Heferina  to  previous  v;ork  on  the  rotatin^^  outer  cylinder 
to  which  the  following  considerations  are  devoted,  first 
A.  Goromerfelds  in'^estlcation  of  couette  flew  between  plane 
parallel  walls  ^  4  j  has  to  be  mentioned.  Although  in  the 
following  Couette  flow  between  concentric  cylinders  will 
be  regarded,  tills  work  has  to  be  mentioned  as  pleine  parallel 
walls  represent  the  limiting  case  of  a  small  gap  between 
concentric  cylinders,  Sommerfeld  proceeded  to  the  equation 
which  determines  the  eigenvalues  but  did  not  analyse  or 


even  discuss  this  equation.  The  only  comment  given  to  this 
equation  Is  that  It  yields  the  time  constant  which  Itself 
Is  expected  to  determine  critical  Reynolds  numbers.  Sexl  [5  ] 
gave  the  analogue  derivation  for  centric  cylinders.  The 
equation  which  yields  the  time  constant  is  the  same  as 
derived  in  the  following.  But  also  here  no  solution  of 
this  equation  was  attempted  obviously  because  of  the  com¬ 
plicated  functions  occuring  in  this  equation  which  are 
similar  to  Sommerfelds  case  Bessel  functions  of  complex 
order  and  argument.  L.  Hopf's  investigation  on  plane  Couette 
Clow  [  6  j  may  be  left  aside  as  his  boundary  condition  of 
coiastant  pressui'e  meens  a  free  surface  and  not  a  rigid  sur¬ 
face  as  it  is  presented  by  the  cylinders  of  the  Couette 
apparatus.  Finally  von  Mtseq  treated  the  stability  problem 
of  the  plane  Couette  flow  as  an  oscillation  problem  [7]  • 

The  eigenvalues  of  a  pai'ameter  are  determined  by  the  zeros 
of  a  polynomial.  Introducing  approximations  only  real  values 
of  this  pararaetei’  were  found  whilst  the  physical  meaning 
of  it  demands  an  iniaglnary  quantity.  This  result  was  Inter- 
pretated  in  the  ivay  that  Couette  flow  should  be  stable  for 
all  nodes  and  Reynolds  numbers  [  8  ]  .  I.ut  this  is  in  con¬ 
tradiction  to  Coaettes  experiments  [^9]  as  well  as  to  raore 
recent  experiments  performed  by  Mallock  [10]  ,  Taylor  |^1l  |  , 
R'endt  [12J  .  Su'jimorfold  regards  this  contradiction  as 
l  -iportant  enougli  to  ask  In  his  book  [15]  if  the  well  established 
iflcchod  of  small  oscillations  should  uDt  be  applicable  or 
if  finite  disturbances  have  to  be  introduced  or  if  even 
the  Navler  Stokes  equations  would  be  insufficient.  Later  a 
more  cautious  interpretation  of  the  theoretical  results 
was  given  by  Lin  [  ]  which  expresses  that,  a  complete  stabi¬ 

lity  proof  still  is  missing. 

Recent  oxnerlnental  results  [15]  hov.»ever  showed  that  the 
earlier  observations  should  be  regarded  with  scepticism. 

Indeed  it  vias  found  that  turbulence  at  least  not  needs 
to  occur-  Objections  on  the  method  of  obsor'/atlon  [  I6  ]  of 
thJ.s  result  may  be  put  aside.  It  was  argued  that  the  inuaersed 


;> 


purtlclOG  which  made  the  flow  visible  should  have  been  sedl- 
mentated  by  centrifu^ai  forces  so  that  no  turbulence  could 
be  observed.  But  attention  had  been  payed  to  tills  phonoraenon. 
bcdlraontation  clearly  could  be  observed  and  turbulence  if 
it  had  occui’od  v;ould  huvo  prevented  sedimentation. 

The  result  that  Couette  flow  with  rotating  outer  cylinder 
may  bo  stable  lead  the  attention  to  side  effects  that  could 
have  Influenced  the  earlier  experimental  results.  Indeed 
it  could  be  shown  experimentally  [l5j  arid  theoretically  [17] 
that  hidden  exccatrlcltles  eind  vibrations  con  be  responsible 
for  the  transition  to  turbulence.  The  theoretical  invest!- 
i^atlon  showed  that  this  Is  not  an  Instability  but  a  seperu- 
tlon  effect  occui*lng  beyond  certain  Keynolds  numbers,  l-iore 
recently  in  an  already  mentioned  Gottinger'  thesis  [thj  atten¬ 
tion  was  payed  to  end  effects.  A  heavier  liquid  was  in¬ 
serted  as  to  seperate  the  teat  liquid  from  the  lower  end 
of  the  apparatus,  obviously  the  feature  of  a  ball  bearing 
was  attributed  to  the  heavier  liquid,  nut  as  slip  neither 
occurs  at  walls  nor  at  contacting  surfaces  of  liquids  by 
this  provision  in  no  ways  the  end  collect  is  annihilated 
but  merely  reversed^ 

This  rcvlev;  ■.'n  recent  'vork  shows  th'.t  the  intrinsic 
ctabllity  problem  of  Couetto  flow  \Jlth  the  view  to  the 
transition  to  turbulence  Jet  is  not  solvou  i.atlsfactorily . 
There  remains  to  bo  deduced  a  conclusive  theoretical  stabi- 
lit:/  '^roof  excluding;  any  presumptive  side  oifccts  on  the 
basis  of  a  Faylelgli  perturbation.  This  will  be  done  in  the 
follov/ing  in  the  f.uamo  work  of  a  correct  ilnoai"*  pertur¬ 
bation  theory.  The  equation  which  yields  thv  eigenvalues 
is  as  was  mentioned  before  identic  with  the  one  derived 
by  Bexl  [  5  ]  •  This  equation  will  be  analysed  by  talcing  full 
account  of  the  functions  oocurlng  In  it. 

Pefering  foi*  the  sake  of  cojaoletonoss  to  the  second  problem 
characterized  by  a  rotatlag  inner  cylinder  and  the  outer  cy¬ 
linder  at  rest  Taylors  calculations  showed  excellent  a':'ree- 
raent  with  observations  [  2  j  . 


As  waa  mentioned  before  two  kinds  of  perturbations  were 
Introduced  In  the  stability  theory  of  couette  flow.  To  begin 
with  Taylors  rotational  syrmnetric  cellular  vOrtwA  pertur’’ 
batlon  the  following  notations  will  be  used:  y  the  radial 
wall  distance;  z  the  axial  coordinate;  u«  v,  w  the  tangential 
radial  and  axial  perturbation  velocity  components;  t  time 
and  fb$  6  real  constants.  Then  this  perturbation  Is  expressed 
by 

u  =  u.  (y)  cosy  z 

(it  ^  * 

V  a  e'  (y)  cosy  z 

w  =•  w^  (y)  slny  z 

It  readily  may  be  seen  that  the  continuity  equation  Is 
satisfied.  P 0  iioans  stability  andP>0  instability.  On  this 
perturbation  bases  G.I.  Taylors  stability  criterion  for  a 
Couette  flow  with  I'Otating  Inner  cylinder  and  vjsi  outer  cylinder 
at  rest.  Application  to  the  reverse  case  of  an  Inner  cylinder 
at  rest  and  an  outer  cylinder  rotating  did  not  yield  results 
[  2  ]  so  that  for  the  present  purpose  this  perturbation  Is 
not  the  proper  one. 

The  other  kind  of  perturbation  was  Introduced  by  Ray- 
lel^  [  18  ]  .  For  Its  representation  Instead  of  the  polar 
coordinates  r,  there  will  be  used  the  distance  y  from 
the  outer  wall  with  positive  direction  to  the  center  and 
the  circumferential  coordinate  x  refering  to  the  outer  Wall 

X  =-  ,  (1 

where  the  radius  r^  is  related  to  the  outer  cylinder*  Then 
the  pertvirbatlon  is  represented  by  the  atreamfunotion 


»  (p  (y)  e 


Y 


i(«AjC  -  /it) 


(2 


is  u  retvl  and  p>  a  nonplex  constant i 


</^  la  related  to  the  wave  length  X  by 

^  -  f 

and  to  the  number  k  of  waves  on  the  olrcunfercnce 


k  -<^r^  (5 

/I  p  Is  the  angular  frequency  of  the  perturbation  whilst 
means  stability,  Instability. 


The  complex  velocity  of  wave  propagation  is 


0  =  Gj.  i  Cj. 


cA  ^ 


(6 


It  la  related  to  r^  according  to'  the  moaning  of  x.  The 
angular  velocity  of  these  perturbation  waves  with  radial 
fronts  la 


This  perturbation  will  be  Inti'oduced.  It  may  be  mentioned 
that  la  this  way  recently  the  stability  of  a  liquid  ro-* 
bating  as  rigid  body  could  be  solved  in  closed  form  in¬ 
cluding  all  modes  [  19  ]  .  The  eigenfianctions  were  found 
to  be  the  Bessel  functions  of  the  first  kind*  As  the  zeros 
of  these  functions  are  located  on  the  real  axis  there 
exists  stability. 


-  6  - 


3. )  Basic  Equations 


Introducing  (2)  in  the  Navler  Stokes  equations  by 
taking  full  account  of  the  viscous  terms  and  limiting 
to  linear  inertia  terms  one  obtains  the  differential 
equation  [15] 


r,  2  5 


r.-y 


r  -y  U'  U 

(u  -  c)  ♦  B  (-U"  +  7^  +  Xr;=y7 


2) 


lA. 

d. 


V.  2  u  2r*  o 

(7^) 


t  t 


2r. 


(ri-y)' 


cp’- 


4r 


(Ti-y) 


,11 


9'  _  <P"  2 

r^(r^-yj^  ’  r|(r^-y)  ~  ^ 


t  t  I 


(8 


•0  Is  the  kinematic  vlaoosity.  U  denotes  the  velocity  of  the 
uoUisturbed  Couetto  flow.  If  V  Is  the  velocity  of  the  crater 
cylinder,  the  .radius  of  the  iVLJo.r  cylliider  one  has 


If 


1 


rv  2 


r.-y. 


1  : 

r^-y 


(9 


Tb'UC  the  velocity  field  is  a  supernpsitlon  vbf  a  rigid  t^djr 
rotation  and  a  hoteritlal  vortex.  Introducird^*  (9)  In  (8); 
and  making  use  of  the  operator 


L(y)  »  ’ 


9 

r^-y 


-e#. 


r.  2 

(77V)  - 


( I0 


Oiie  obtains 


fji  - 


1  0 
“3r 


r,-y 


L) 


(11 


VJith  (10,11)  the  fourth  order  equation  (8)  is  reduced  to 
tvio  equations  of  second  order. 


Introduelng  dlinenalonlaae  quantities 


5. 

g 


Vr. 


R 


and  t-r'isnsf oy^  t'. 
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(  1“  1j  )*"  f  ’  ’  “  (  1"-  1|  )(?’-  k* 


r’‘:j  { 1“  ^  )"L 
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7-  '’-iJM 


i 

FT? 


(1-  ij)!- 


*  1 


L'  - 


Now  Introduolnc  (9)  and  the  abbreviations 

1 


B  - 


■f:'? 

(yl)  -1 


3  P 

-  1  k  R  B  ^ 

r  2 

b^-  ikRt-llcRB(^)  , 

o 


(  12 


(1? 

(H 


(15 


(16 

(17 


further  Introducing  the  aubatltutlon 

a  -  b  d--^  )  ,  (18 

(14)  tranaferM  If  multiplied  by  k  R/l  to  the  Beaael 
equation  [ 13  ] 

a^L’‘  +  a  L’  +  (z^-5>^)L  -  0  ( 19 

Intreduoing  the  oylinder  f\mctlon  of  oonplex  order  q 
the  aelutlon  of  (19)  ean  be  written 


Zj,  (*) 


(So 


Tti«@rtSng  this  solution  1?a  ( I3)  one  obtains 


1 1 


(z) 


The  Bolutlon  of  The  howogenous  part  of  (21)  Is 


(?1 


<p^  -  z  z  .  (22 

The  particular  solution  of  (21)  la  feunu  by  th#  ■othod 
of  variation  of  constants.  Satisfying  the  bounosry  con¬ 
dition 

Tp  “  ■  P  Oq)  "  0  (23 

one  has 


z 

r 


-k-k  _k  c-k 

z  4  -  z  4 


(24 


Here  correspO'i'Is  to  :■  0  what  meCLii3  according  to  (18) 


2o  b 


(25 


Prom  (24)  follows 


Tp  “  -rr  (5)  <14 


(26 


With  tho  notations 


»  z 

,(z)  -  /  4“*'  Z,  (4)  d4  ,  I./Z)  -  J  4*'‘*’  Z  (4)  i  4 


(27 


the  partioulsr  solution  is 


fp  -  - 


1 

7F 


-k  k 


^6  <l#}f’l¥atl'f'Qs  tilth  raepeat  to  s  of  the  integrals  in  (26) 
will  be  Oao  baa 


d 

n 


[  Z^(C)dC 


Z  (s) 

P 

-1 


d 

*  dr 


;  f  s""*’  (5)a  t 


■‘'*’z„(*) 


/ 


tilth  tbfiSfi  oApfoaiSlOila  Gilo  uGFIToS  rpCS  (c6)  )dl6Q  Intl^* 

duolng  (27) 


♦p  -  b'*  *“  b)  • 


The  eOBplAte  aolutlon  and  its  derivative  are 
,  .  C,  z*'  .  Cj  z-''-  ^  I,  -  Ij) 

k  /-  k  -  „"kv  .  1  .-k~1  T  .1  _k-1  T 

9’  -  -  (C^  z  -CgZ  )+^Z  j  I^+J*  Ig 


lntz*oduolng  the  boundary  oondltlons  9  «  f'  0  at 
one  has  If  tntroduoing  cpp  ■  0  at  thia  boundary 


*&*'  *  ®2  *0  ■  ® 


®1  -  =2  =0'“  -  0  • 


vhat  neana 

*  Cg  -  0  .  -  <28 

Further  the  same  boundary  oondltlons  at  yiold  when 
regarding  (26) 

*•1’^  *1^  ^2  <®l^  “  ° 

a,’*'  I^(z^)  ♦  z^  (z^)  -  0 


«  “jQ  ^ 


sqimt-iofis  ■©an 


■k; 


>  1  ^  ^  V 


V  ^  ^  u  fi 


^  ^  i  \ 

>  “I  *^2  ^  ^  ^  ^ 


\j 


hm  is  &n  arbitral^  b-otsidaiTr  ther*e  euiet  be 


I,(ap  -  0  .  .  I2^*1^  “  ° 


(39 


Nom  aepar&tiois  the  o^llnder  funotlon  In  Ita  two  llnvairllj 
indepwident  parts  N*H  with  oooatants  C^« 


Z„  -  C,  N  4-  Cj.  » 
f  3  ^ 


oue  obtains  froa  (37«  29) 

Z. 


I’ 


Cj  1  N(4)  d5 


M(J;)  d  ^  -  0 


^1 

I  N(t)  I  M(c)  dc 


As  are  arbitrary  the  determinant  axist  be  sero.  This 

condition  gives 


1  1  1 
j  5“*’  N  4  5  •  f  M(t)  d  t  -  J 


M  d4 


i’ 


»  dc 
<3o 


whioh  Is  the  relation  derived  by  Sexl  for  olrowoferantlal 


Couette  flow 
Couette  flow 


5 

K 


and  before  by  Scaaerfeld  for  r»otlll3iear 
•  Ihe  boundaries  are  aooordlng  to  (18,  25) 


If  la  written  for  ^ 


*o  -  to  »  s^  -  b  (1-  tj) 


(31 


“  11  - 


lo  that  (3^)  should  be  satisfied  at  least  by  one  pain  of 
falues  bg 'f|  B  It  Is  seen  that  (50)  is  aatlefled  for  arbl=- 
tram’  'Walues  n  b  ®  0.  uould  ^an  neutral  atabi  - 

I  ^  -  * 

lity  and  a  wave  propagation  with  the  rigid  body  ^rt  of  (9). 
But  ^o,  ~di)  show  that  b  «  0  would  not  yield  a  perturbation 
aa  f  would  be  independent  frow  L  wtiat  oont radio  ts  to  the 
Meaning  of  L.  Thus  b  »  0  has  to  be  excluded. 

^.)  Tranaforaatloa  of  the  tranaoendenfcal  equation  of 

80— erf«ld~3exl 

To  analyse  (30)  series  expansions  of  the  integrands 
starting  at  will  be  Introduoed.  Writing  H  for  M,  N  one 
has  with  (31) 

tk+1  +k+1  -flc+l  , 

4  H(^)  -  b”  H(b)  +  (r  H(t)  )  (C-b) 

^-b 


+  b(^)  )”  —  (&-b)2  I-  . 

^«b  2! 

By  integration  one  obtains  with  (31) 

p  H(t)  dt  -  H(b)  +  H)  .  (32 

z  ^ '  5-b 

o 

'Equating  the  ooefflolents  one  obtains  when  writing  k  for 
±  k  and  H  for  H(b)  as  to  slspllfy  the  notations 

H)'  -  (k+1)  b*^  H  +  H' 

4-b 

jjj"  -  (k+1)  k  b*^“^  H  +  2(k-»>1)  H*  b**'*’^  H*  ' 

C-b 

H)’”  -  (k+1)k(k-1)b*'*^  HO{k+1)k  b^"**  HV3(k4l)b‘^' V 


1^?  ™ 


Heax’ranglng  noif  the  terea  of  (.T^)  with  rsapeot  to  the 
derlvatlvaa  of  H  one  obtains 


^1 
f  * 


j  H  ri,4  =  H*'  +  •••. 


U 


The  coeffloients  are  deteralned  by 


(33 


oo 


n- 1  k+ 1 


)(  ) 

n!  n-^-1 


n*  V)  +1 


(34  a 


Introducing  (33)  iii  (3®)  and  the  notation 


a(lc)  =  a  ,  a(“k)  c 


(34b 


one  obtains 


(Md  ' -M' M)  (a^c  )+b^(MiJ '  ' -M’ '  T)  (a^Cg-agO^) 

+  b'^{(MN' "-M*  ' ’>0(a^c^-a^c^)+(M'N' )(a^C2-a20^)} 

+b®{(MN^^-M^\')  (a^C;^-aj^OQ)-»-(M’N'  "-M'  *  'n'  )(a^o^-a^c^)  J  (3^ 

+  (M"N»  "-M»  ”  r  • )  (agO^-a^Og)]  \ 

+b‘'°{  (MM'^-J^Nl)(aQ0g-agCQ)+(M’N^-M^N')(a^05-a50^) 

+  (M' » )  (ago^-a^og)! 

+ . -  0 

After  dividing  by  this  equation  is  represented  by 


l>«0 


r+1 


(!“  (n)  )_„(r“n+1  )j.(b)  w 


a  o  ^  “ft  jQ  )  s>  0 

n  r-n+1  r-m-l  n 


n«0 


(35b 


whloh  is  the  transformed  Sonmerfeld-Sexl  equation  (3o)* 
This  new  oxpresalon  permits  as  will  be  shown  the  direct 
numerical  ewaluation  of  the  stability  parameters 
for  all  modes  k  with  any  desired  exactness i 


5 • )  Solution  of  the  fundamental  equation  (35b) 


The  Wronski  determinant  W  will  be  Introduced.  One 

he.s 

w  -  MN'  -  NM’  ,  w”  -  nr;”  -  nm” 

With  the  abbreviations 


z  z 

the  equation 

y”  +  f  y'  +  g  »  0 

has  the  solutions  M,  N  according  to  (19). 
■nie  differential  equation  for  W  is 

W 

- - 
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(36 


(37 


(38 


(39 


and  therefore 


-  u  - 
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?  $ 


r’u  a.  «  ---  W 
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f  t  I 


W(-  f • '  +  ^  ff ’  -  f 


-  r^\ 
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»  W(-  f’"  ■»■  3  f'^+  4  ff’'~6  f^f'  +  f^) 


24 

7“ 


(4o 


Vf’ 


W(-  +  1o  f'f”  +  5  fl'”'  -  ir  -  10  f^f" 


+  1o  Pf'  -  -  - 


12o 


Prom  this  one  deduces 

=  (-1)®  (4l 

2® 

On  the  other  hand  It  follows  from  (36,  39) 

;  -  v;”  -  (M’N’‘-M’'N’) 

-  W’"  -  2 

(42 

MN^  -  m'^N  -  -  2  (M"N" ’-M" 'N”  )'-3  ( M * '  ) 

-  2  (M"N^^  -  M^'^N’')-3  ' ) 

-4  (M*N^-m'^N*) 

Inserting  fJ,  M  in  (38)  and  fonaing  the  m-th  derivative 

m(«)  +  (f  '  *  ^  g  -  0 

+  (f  N'  +  g  -  0  , 


and  multiplying  the  first  equation  with  the  n-th  derivative 
of  Kf  and  the  second  with  the  same  derivative  of  N  one 


obtains  whsn  subtraotlng 


(f  M'4^g  M) 


(f  N’+g  N)'“  -'n'"'  j 

/  h  •« 

v-t^ 


with  thle  expression  the  higher  derivatives  oan  be  ex** 
pressed  by  the  lower  ones.  Finally  the  Wronskl  detenalnant  W 
can  be  Inti^duoed  so  that  In  the  baslo  equation  ov  33b) 

all  functions  can  be  expressed  by  W.  Thus  In  (35)  all 
functions  M,  N  and  their  derivatives  are  replaced  by  W, 
which  shows  to  be  a  ootsaon  factor  and  therefore  oan  be 
canceled  In  (33a)  und  (33b)  reap. 


Indeed  fron  (43)  one  deduoes 


N*  ’M  -  f4’  'N  »  -  f  W 
N**M*  -  M"N‘  »  g  W 

-  M”N”*  -  (fg'-gf’-g^)  W 
-  M^^M”  »  (ff ’c+B^f“f^g+fg"-Gf  ’  »-2gg')W 


N'''H  -  M”*?!  (-f’-g+f^)  W 

N”'M'  -  K*"N'  •  (-f»^g')  W 

-  M' **?:’*  -  (f ’g+g^-fg’ )  W 

-  (“f  * '-agV^ff '+2gf-f^)  W 
-  M^^N'  -  (-af’g-g^-ff^g-fg'-fc’ ‘  )  W 


(36,42)  Introduced  In  (35a)  gives 
b°  tf(a^o,-a,o„) 

+b’  Vf'tajjOg-agO^) 

+b®  jw’  •  (a^Oj-ajO^)-(M'»"-M'  ’N' ) 


.>1) 


-  16  “ 


(y ' *  ^ ' ’-M”  'H' )  [  2(a^o^-a^0^)-(a^o^-a^3^ ) 

{w^^(a^Oe-ftKO,J-(M’  'n"  ’-M*  '  'N* '  )f  2(a^0g-a^o„)-(&^o,-a,0g) 


r.  -a  rt  "Ja^a  #!.-«./»  \  11 

...  ..  ..  '-1-4  -4-rjj 


+b' 


5(a^og-ag0Q)-{a20|^-a^O2)j 

] 


- (M*N^-M^N * )  "5"o'  '"1''5  “5"1  j 

{w'^(a^o^-ayOQ)-(M*  ’  ’ ' )  j■  5(a^o^-a^o^)-(4^04■-^^«^^ 

-Ym* ‘ 


(44)  Introduced  in  (45)  gives  with  (4o)  if  the  oo««on  ■ 
factor  W  la  cariceled 

+b^{(a^,o^-ajCo)(2f2)»  [  (aoC^-a^o^l-Ca^Og-agC, )  ]  g  j 

+b^{(aoO^-a^Oo)(-6p)-  [2(ao04-a|0^)-(a^C5-a  o,)l  (-f*fg')} 

4/  k  '  '  X46 

+b  l(a^05-a50^)(24  f^)-  |  2(a^,03-a^c^j)-(a2C2-aj02)  j  (f'^^-fg*) 

■  5(a^05-a^o^)-(a^C4*a^o^)]  (-2f 's-g^+f^gl 

J  -fgV^g**)  J 


4 


0 


•  17 


In  f  and  g  the  boundary  ^  b  haa  to  be  inserted  for 
7,  aooordi.ng  to  {'^^)  *  Thla  meana  that  in  the  first  ooluron 
of  (46)  b  oanoels.  Further  It  has  to  be  raentioned  that 
aooordlng  to  (34a)  the  first  brackets  In  the  two  oolumna 
ui  (46)  iiioxeaae  with  pOwei'a  uf  n  •  Tlie  lii'at  uOlUuiii 
begins  with  the  first  power  and  the  second  with  the  third 
power.  Thus  (46)  represents  a  series  expansion  in  powers 
of  Tj  , 

6. )  Dlsousalon  of  the  functions  aj  (n)#  o.)  (n) 


One  haa  according  to  34 (a, b) 


oo 


n«-0+ 1 


n 


nl  n->J-1 


+k+ 1 

)(-  ) 


(47 


This  expression  shows  that  the  functions  are  polynomials 
and  that  c^  are  series  with  terms  negative  for  even  in¬ 
dices  and  positive  terme  for  the  odd  ones.  One  finds 

that  are  the  futict.lons 

o  o 


1  1  +k+2 

- - +  — —(1-4)  (48 

-tk+S  -tk+2 

This  expression  shows  that  is  negative.  For  c^  one 
has  for  kic.1 


°o  “  ■  ^ 

•nie  undetermined  expression  (48)  for  k  **  -2  is  equal 
to 


o 


o 


In  (1-T|) 


ia  " 


For  one  dcriveo  from  (^7) 


1 

1-t] 


The  a' a  show  to  have  first  ascending  and  then  descending 
terms-  For  odd  Indices  the  sum  of  the  ascending  terms 
la  smaller  than  the  largest  positive  terra  and  for  even 
indices  this  sum  is  smaller  than  the  largest  negative 
term  If  In  each  case  the  largest  terras  are  excluded 
from  the  addition.  As  the  descending  terras  are  descending 
strongly  and  as  there  number  Is  limited  the  a' a  with 
the  odd  indices  show  to  be  positive  and  with  the  even 
Indices  negative.  This  statement  holds  for  all  0 «  *  1 

so  that  In  these  regions  the  a' a  have  no  zeros. 

Table  I  gives  the  values  of  a,  c  for  k  >=  1,2,3  and 
table  II  the  values  of  the  brackets  ocouring  in  (46). 

All  these  values  refer  to  tj  =  0,5. 

7 • )  Evaluation  of  the  stability 


Inserting  (37)  and  z  ==  =  b  in  (46) 


one  obtains 


-(aoOg-agOo) 

+2(a^c^-a^c^) 

-6(aoC4<-a^Co) 

+24(a^C5-a5C^) 


(a^c^-a3C^)-(a^C2~a2C^)  J 

2  ( aQO^“a^c^ )  “  ( ^  ^ 

- 

3(a^a^-a^c^)-(a^c^-aj^c^) 


(49 

(b2-s2) 

(-b^+3j^) 


0 


-  iy 


Inserting  with  table  II  the  values  for  k  ^  1  the  real 
und  imaginary  part  of  (49)  yield  two  equations, 

p  1  t.h ^  ij«  t:  1  OH 


18,68760851  +  7»6i5938o2  r  7»42o258o2f^ 

-  0,15278165  (b^  0 


{50 


1?  obtained.  With  (15^  l6,  17)  and  the  numerical  value 


^  ^ 


K*"  ,  corresponds  to  the  assumed  value  »  o,5  the  se¬ 

paration  of  the  real  and  imaginary  part  of  (5o)  gives 
the  two  equations 

-  0,15278165  -  7,9215013  R4^  +  11,11458885 

+  0,15278165  -1)^  =  0  , 


(51 


7.9215013  0.0652333 

0,3055633  {^p-D 


(52 


Inserting  (52)  in  (51)  one  obtains  exactly  a  biquadratic 
equation  In  tj,“i  , 

R^(4r-1)^  +  2088. 949384(6^.-1)^  -  0,0455760  -  0 
in  which  R  appears  as  parameter.  The  numerical  solution  is 


R  - 

1 

“  ± 

4,69  • 

1 

0 

R  - 

I0 

6j,-"i 

"  ± 

4,67  • 

1 

0 

R  « 

loo 

«  + 

4,67  • 

R  - 

00 

^  0 

n  -- 


Only  the  nor.atlve  values  have  physical  meaning..  Intro" 
<iMcln(.;  them  in  (52)  one  obtains 


R  = 

1 

-  71,44 

ntr 

--  •-r -1 

11 

%\J 

~  1  <  f 

R  « 

loo 

h4^ 

«  -  71,63 

This  result  means  stability  and  it  la  seen  that  4^  la 
Inversely  proporciorval  to  the  Heynolda  number  R.  The  nunerloal 
values  obtained  for  4p"1  ^ow  that  the  circumferential 
propa;iatloii  velocity  of  the  perturbation  is  extremely  near 
to  the  circumferential  velocity  of  the  outer  cylinder  with 
a  nei:,ll(^able  influence  of  the  Reynolds  number.  This  result 
la  analo,j;ue  to  the  result  found  for  the  stability  of  a 
liquid  rotating  au  rli_,ld  body  for  which  in  our  notation 
r^/r^  =  oo ,  where  4^*1  *=  0  was  foujiu  [  j  • 

For  k  =  3#  T  =  o#5  (^9)  yields  by  neiilectlng  the  last 
terra  of  (46) 

^  o  2 

“  o,o1oyo6  =  -  o, 008207  b  +  o,ooS45y 

r*.eparatin:j  the  real  and  Imaginary  parts  with  (16,  17) 
one  obtains  the  two  squat ions 

-  0,010766  »  0,008207  icR4j^  0,00845 

* 

C,  e. 

-  o, 008207  kR  (4j.  -  P'(~)  -  0,008716  kRB  -  0 

o 

The  first  equation  gives 

R4^  -  -  3,53 

Olid  the  second  equation,  when  introducing  with  (15)  r^/r^  “  2  , 

4j,  =  o,98l  . 


I 


2l  - 


Again  stability  and  a  propagation  velocity  near  to  the 
eireumf erentlal  velocity  of  tha  outer  cylinder  is  found. 

Here  appears  to  be  Independent  of  the  Reynolds  number  R. 
Indeed  the  influence  of  R  appears  only  in  the  next  hij^er 
approximation  that  means,  when  the  last  term  of  (46) 
is  oorisidered.  Therefore  in  any  case  the  Influence  of 
R  on  ^  Is  weak  as  It  was  found  numerloally  In  the  case 
k  ■  1  by  considering  also  the  last  term  of  (46). 

®«)  Coaolualons 

The  Soouierfeld-Saxl  transcendental  equation  which 
yielda  the  stability  parameters  of  Couette  flow  was  first 
transformed  (sue  aq.  33b).  Then  by  Introducing  the  Wronski 
determinant  it  was  possible  to  eliminate  the  unknown  Bessel 
functions  of  iuiaglnary  argument  and  imaginary  order#  whlob 
until  then  prohibited  the  solution  of  the  original  Suiuier'' 
feld-Sexl  equation.  A  series  expansion  in  the  wall  distance 
with  the  stability  parameters  In  polynomial  fora  was  ob*^. 
talned  (see  eq.  49)  which  can  be  solved  numerically  for 
any  modes  and  widths  of  the  gap.  Pur  a  ratio  of  the  Inner 
and  outer  radius  of  2  stability  was  found  for  the  first 
and  third  mode. 


'i'anie  I 


Coeffl 

clents  a;  ,  for^j  «> 

0,5 

% 

k  «  1 

k  2 

k  -  3 

■  h 

c 

_  2 

-  La  2 

-  1 

0 

2 

11 

lo  12 

19 

TOT 

c 

1 

1 

o 

3^ 

lO  -  12o 

563 

24-60 

A. 

-  J- 

.  JL. 

.  -Si. 

1D 

lOo 

1120 

c 

2 

1 

*  w, 

-  113 

330^ 

1447 

35t)- 105 

11 

1 

_^163  ^ 

^3 

o 

CVI 

Iter 

16 -16- 24 -70 

c 

3 

1 

w 

96  ^0 

3071 

29 

31  _ 

1 

96-96-35 

32-32-16-35 

54-16-35 

c 

1 

6887 

1439 

32-120 

32-288-48-35 

32  27-64-35 

Table  II  Coefflclenta  of  equation  (49)  for 
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